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Abstract 

This  is  the  final  report  for  Phase  I  of  the  above  STTR  contract  with  AFOSR  during  the 
period  September  1,  2002  -  September  30,  2003.  It  is  based  on  Status  Reports  1,  2,  and  3. 
We  summarize  the  development  of  integral  equations  and  their  numerical  implementations 
in  solving  wave  propagation  in  a  cylindrical  fiber  with  concentric  dielectric  layers.  These 
multilayered  fibers  are  also  called  Bragg  fibers  [12,  13,  15].  From  preliminary  computational 
results,  we  find  that  an  air-core  multilayered  fiber  can  effectively  guide  wave  in  its  air  core 
when  its  total  radius  is  a  quarter  of  the  operating  wavelength. 


1  Introduction 

An  omnidirectional  reflector  which  reflects  electromagnetic  radiation  from  all  incoming  angles  and 
over  a  broad  range  of  frequencies  but  exhibits  virtually  no  absorption  losses  has  been  fabricated 
for  the  20,000  GHz  -  30,000  GHz  frequency  range  [4].  Recently,  an  omnidirectional  dielectric 
mirror  fiber  which  is  an  omnidirectional  reflector  of  cylindrical  shape  with  the  diameter  of  200 
\im,  has  been  developed  by  Yoel  Fink  and  his  colleagues  [7].  These  omnidirectional  dielectric 
mirrors  can  be  used  in  applications  requiring  optimal  confinement  or  reflection  of  light  at  all 
angles,  such  as  optical  cavities  or  hollow  waveguide.  It  is  of  great  interest  to  the  Air  Force  if 
omnidirectional  reflectors  and  fibers  which  operate  in  the  lower  frequency  range  can  be  fabricated. 
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Microwave  (MW,  1  GHz  -  18  GHz)  and  millimeter  wave  (MMW,  18  GHz  -  100  GHz)  frequencies 
are  particularly  desirable  (see  Fig.  1)  in  radar  technology. 


Figure  1:  Electromagnetic  spectrum 

Based  on  experience  and  published  data,  it  is  likely  that  high  performance  omnidirectional 
reflectors  for  lower  frequencies  (MW  and  MMW)  will  be  too  thick  to  be  practical.  Hence,  at  lower 
frequencies  there  is  a  need  for  higher  contrast  ratio  layer  structures,  to  enhance  the  reflection 
coefficient  at  each  interface  to  maintain  a  reasonable  overall  thickness.  However,  the  use  of  higher 
contrast  systems  will  lead  to  a  reduction  in  angle  coverage.  Therefore,  to  determine  the  low 
frequency  region  where  an  omnidirectional  reflector  and  a  low-loss  dielectric  fiber  are  viable  and 
candidate  materials  for  achieving  required  performance,  it  is  important  to  perform  analytical  and 
numerical  study  of  the  wave  propagation  in  the  planar  and  cylindrical  dielectric  layered  structures. 

In  our  preliminary  study  of  Bragg  fibers,  it  is  possible  to  concentrate  electromagnetic  wave 
in  an  air  core  Bragg  fiber  whose  radius  is  comparable  to  wavelengths  in  the  microwave  region. 
Further  investigation  and  analysis  are  necessary  to  reach  a  more  definite  conclusion. 

Omnidirectional  reflectivity,  the  ability  to  reflect  waves  at  arbitrary  angle  of  incidence,  is 
possible  in  multilayered  dielectric  media  whose  band  diagrams  consist  of  spectral  regions  where 
no  light  transmission  exists  in  the  media,  which  are  called  band  gaps.  For  planar  multilayered 
stacks,  transfer  matrix  method  has  been  used  in  computing  their  band  structures  and  reflectance 
spectra  (see,  e.g.,  [1,  2,  5,  14]). 

Unlike  conventional  fibers,  which  guide  light  through  total  internal  reflection,  Bragg  fibers 
achieve  light  confinement  through  Bragg  reflection  [13].  Bragg  fibers  uses  a  one-dimensional 
periodicity  of  concentric  rings  that  exhibit  omnidirectional  reflection  in  the  planar  limit.  In  other 
words,  in  the  limit  of  large  radius,  the  circular  Bragg  fiber  approaches  a  planar  Bragg  stack,  whose 
band  diagrams  consist  of  band  gaps.  Therefore,  it  is  necessary  to  have  Bragg  fiber  modes  to  fall 
inside  the  TM  and/or  TE  band  gaps  of  the  corresponding  Bragg  stack. 

Wave  propagation  in  a  Bragg  fiber  can  be  analyzed  by  a  transfer  matrix  approach  and  the 
minimization  of  radiation  loss  as  in  [15],  or  by  numerically  solving  Maxwell’s  equations  using 
finite  difference  method  in  both  frequency  and  time  domains.  Recently,  asymptotic  matrix  theory 
of  Bragg  fibers  is  developed  in  [13].  Yeh’s  transfer  matrix-minimization  method  is  complicated 
for  guided  modes  that  are  not  TM  nor  TE.  In  order  to  use  finite  difference  methods,  one  needs 
to  construct  fictitious  boundary  conditions  such  as  perfectly  matched  layers  (PML)  or  absorbing 
boundary  conditions  (ABC).  They  can  be  also  time-consuming.  The  asymptotic  matrix  analysis 
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fails  if  the  fiber  core  becomes  too  small  and  it  is  difficult  to  estimate  the  accuracy  of  the  asymptotic 
results.  Operating  wavelengths  in  these  articles  are  in  the  optical  region  of  the  electromagnetic 
spectrum. 

'  The  rest  of  the  report  is  arranged  as  follows.  In  Section  2,  we  present  the  electric  field  integral 
equation  which  describes  wave  propagation  in  a  cylindrical  waveguide  of  arbitrary  shape.  In 
Section  3,  Nystrom  method  is  employed  to  discretize  integral  equations.  The  approach  is  simple 
but  very  efficient.  In  order  to  determine  propagating  modes,  non-linear  eigenvalue  problem  needs 
to  be  solved.  Newton’s  method  and  Ridder’s  algorithm  are  implemented  to  find  solutions  of 
the  eigenvalue  problem.  In  Section  4,  numerical  examples  are  given  to  illustrate  the  validity  of 
the  integral  equation  formulation  and  numerical  techniques.  Finally,  the  report  is  concluded  in 
Section  5.  Here,  we  also  draw  attention  to  mathematical  issues  that  require  detailed  and  rigorous 
investigation. 


2  Integral  equations 


In  this  section,  we  briefly  present  the  EFIE  to  describe  the  surface- wave  propagation  in  a  dielectric 
cylindrical  waveguide.  First,  the  general  3-D  EFIE  for  a  waveguide  with  a  graded-index  core  of 
arbitrary  cross  section  shape  is  introduced.  The  3-D  EFIE  is  then  reduced  to  a  2-D  EFIE  when 
the  index  profile  n  is  axially-uniform,  that  is,  n(x,  y,  z)  =  n(p,  0).  Finally,  for  a  graded-index 
circular  fiber  with  the  radially-graded  index  core,  i.e.,  n(x,y:  z)  =  n(r),  we  further  reduce  the 
2-D  EFIE  to  a  system  of  1-D  integral  equations.  We  are  especially  interested  in  natural  surface 
wave  modes  which  are  nontrivial  modal  solutions  to  the  homogeneous  EFIE  when  the  incident 
field  El  =  0. 

In  integrated  dielectric  waveguide  system  with  a  dielectric  waveguide  whose  refractive  index 
profile  is  described  by  n{x,  t/,  z),  lying  on  the  planar  interface  {y  —  0}  between  a  substrate  having 
index  ns  and  a  cladding  overlay  of  index  nc  (e.g.,  air  with  nc  =  no),  the  fundamental  electric  field 
integral  equation  can  be  expressed  as 

E(r)  -  (I fcc2  +  VV  )  [  5-^p-C (r,r')  ■  E{r')dV'  =  Einc(r), 

JV  nc 

where  5n2(r )  =  n2(r)  —  n2,  V  is  the  volume  where  Sn2  ^  0,  kc  =  ujy/JIoel  =  = 

konc ,  and  Q(r,  r')  is  the  3-D  Hertzian  potential  Green’s  dyadic  defined  as 

G{r,r')  =  &,{ry)  +  Gr{r,r'), 


with  &  is  the  principal  Green’s  dyadic  for  primary  wave  of  Hertzian  potential  excited  by  polar¬ 
ization  and  Gr  is  the  reflected  Green’s  dyadic  for  reflected  wave  of  Hertzian  potential  excited  at 
the  substrate-overlay  interface.  Note  that  if  ns  =  nc  (the  substrate  is  removed),  GT  =  0.  The 
principal  Green’s  dyadic  is  defined  as 


5p(r,  r') 


ikc\r-r’\  = 

- - /. 

47r|r  —  r'\ 


The  formula  of  Gr  is  omitted  for  brevity.  In  the  rest  of  the  report,  we  are  interested  in  natural 
surface-wave  modes  which  are  the  solutions  to  the  homogeneous  EFIE  when  Ernr  =  0  and  for 


ns  —  tic. 
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2.1  EFIE  for  axially-uniform  waveguides 

Consider  the  2-uniform  waveguide  whose  refractive  index  profile  is  ^-independent,  that  is,  n(r)  — 
n(p),  where  p  —  xx  +  yy.  Consequently,  the  electric  field  E  can  be  written  as 

E(r)  =  e{p)e~t/3z,  p  =  xx  +  yy , 


where  (3  is  the  propagating  constant  which  is  an  unknown  eigenvalue  parameter.  The  homogenous 
EFIE  becomes 

e(p)  -  ( k2c  +  WO  /  c/5'^-^Kp,  p')  •  e(p')  =  0,  (1) 

Jn  nc 

where  f2  is  the  cross  section  of  the  waveguide  and 


=  f  ~  _  ^t\T  K0  hc\p  p  I)  J 

g{p,  P  )  =  9{p,  P)I= - 2^ - 


with 

7c  =  \/P2  -  kl-  (2) 

By  using  divergence  theorem,  one  can  show  that  transverse  components  of  (1)  are  independent 
and  hence  yield  a  unique  description  of  the  surface- wave  modes.  Therefore,  it  is  sufficient  to 
consider  the  following  equation  for  the  transverse  components  et  =  exx  +  eyy  of  e: 


et(p)  ~  Vt  £  5n ^2^  '  et(p')ff(p>  p'W  +  jf  —  ~g(p’  p')dS') 


n2(p) 

+kl  [  Sn2{p')et{p')g{p,  p')dS'  = 


(3) 

0,  Vp  €  ft, 


where  V  =  Sf2,  n!  is  the  outward  normal  of  the  boundary  T,  and  V*  = 


2^  for  axially-uniform  circular  waveguides 

Assume  that  the  index  profile  of  the  waveguide  only  varies  radially,  that  is,  n(p)  =  n(r).  In 
(3),  we  note  that  n'  •  et  —  rf  •  et  =  er,  thus,  the  homogeneous  EFIE  in  polar  coordinates  can  be 
expressed  as 


,  .  /  i  Sn2(r')  ,  ,  Koiidp  -  p'\)  ...  f  dn2(r')  er(r',  </>')  Kofrdp  -  p?\)  , 

e,(r,  -  V,  (-  jl  -^-er(r  ,  *  ) - ^ - d‘  +  Ja  ~d ~ “rfVj - £ - iS 

Sn2(r')e,(r',  — ■«' 


=  0. 


We  can  decompose  the  above  vector  EFIE  into  two  coupled  EFIE’s  for  the  scalar  fields  er  and 
as 


,  ^  d\  1  8n2(a)  ,  j.^K0(jc\p- p'\)  JJt  t  f  dn2{r')  er(r' ,<£')  ICo(7c|p  -  p'\)  / 

-fco  f  Sn2 fr')[c,(r', 4>') cos{<p-  <p’)  +  e.;,(rf  ,■/)  sinf./.  -  ^')]  — — dS'  =  0 

•/  Q  ; 
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«*■*>  - \h  H 

-kl  f  5n2(r')[-er(r',  4>  )  sin (<p  -  <j>')  +  e^(r;,  </>')  cos (<f>  -  — —dS'-  0 

JQ 

We  can  write  the  unknown  scalar  field  components  are  represented  by  Fourier  series  with  radially- 
dependent  coefficients  in  the  angular  variable  <f)  as 


oo 

er(r,<f>)=  Ar ,m{r)el™i>, 

m=— oo 


oo 

iA4>,m{r)eim4> . 

m=— oo 


(4) 


Substituting  (4)  into  the  integral  equations  for  er  and  yields  the  integral  equations  for  the  radial 
coefficients  Ar^m{r)  and  A^m(r):  for  0  <  r  <  a  where  a  is  the  radius  of  the  circular  waveguide, 


AT,Tn{r)  fo  Ar)Tn(r') 


6{P  -  a)6nH/)  1  rfn2(r')~|  dW0i7rn(r, r')^, 


'  ■*' '  [  n£  n2(r')  dr?  J  Ur 

kl$Z5n\r')[Ar,m{v')Wc„m{ry)  -  A^m(r')Wsam(r,r')r'dr'  =  0, 


dr 


and 


A^y-^foArAr') 


8{r'  —  a)5n2(r')  1  dn2(r') 


ni 


n2{r')  dr' 


Wo,7 m(r,r')r'dr' 


-k2  /0“  8n2{r')[-Aryr')Wsnm  (r,  r')  +  A^m(r')Wcr/m(r,  r')}r'dr'  =  0, 
where  7m  is  defined  by  (2)  and  the  kernels  Woi7m ,  Wc>7m  and  Ws>7m  are  defined  as 


W0>7, m(r,r')  =  ^-  J  e^KoAmi/r2  +  r'2  -  2  rr'  cos  ip)  dip 
_  /  Hn(^ymr)In(^ymr  )  if  r  <  v, 

~  \  Im{lmr)Km{lmr')  if  r'  >  r, 


(5) 

(6) 


Wc>7m(r,  r')  =  J  eim^  cos  ipK0  \Jr2  +  r'2  -  2rr'  cos  ip) dip 

jj  4"  Km—l{')/mr)Im—l{^mf'  )  if  T  <  V, 

~  2  \  Im+i  {lmr)Km+ymrr)  +  /m-i  (7™T  )^m-i  (■ymr')  if  r'  >  r, 

and 

W*,'Kn(r»r0  ~  ^  f  eimi’ sin  ^K°  (7m  \A2  +  r'2  -  2rr'  cos  tp)dijs 

_  ~  Km—  li'lmrflm—  lilm^  )  if  ^  ^  **, 

2  \  Im+1  (7mr )-^"m+l (7m ^  )  fro— 1  (7m^)-^m— 1  (’Tm*'  )  if  r  >  f. 


3  Computational  method 

Integral  equations  (5)  and  (6)  are  discretized  by  Nystrom  method  [10].  Since  kernels  Wo,  Wc,  and 
Ws  are  in  closed  forms,  the  implementation  of  the  method  is  simple. 
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3.1  Nystrom’s  method 

Nystrom  method  is  used  to  discretize  (5)  and  (6).  For  brevity,  the  interval  (0,  a)  is  divided 
into  segments  of  length  A r  =  a/N  where  N  is  the  number  of  subintervals.  Denote  r\  =  %  *  h, 
i  —  1, 2, . . . ,  N  and  ro  =  e  for  some  0  <  e  <  1.  First,  the  integrals  in  (5)  and  (6)  are  approximated 
by  the  composite  trapezoidal  rule.  Hence,  we  get 


N 


Ar>m  (r)  ~  WjAr,m (rj  ) 


3=0 

N 


5(rj  -  a)An2(rj)  1  dn2{rj) 

n2  n2(rj)  drf  J 


dW^J^rj)  A 
dr  ^ 


-k20  WjAv?  (rj )  [AriTn  (rj )  Wc,^m  (r,  rj)  j)Ws,ym  (r,  rj)rjAr  —  0, 

j= o 


and 


A(f>im(r)  —  -  Ar>m(rj) 


N 


N 


3=  0 


£(77  —  a)  An2  (77)  1  dn2(rj) 


rtf 


n2(rj)  dr'  J 


An2(rj)[  rj)  +  rj)]rj^r  —  0) 

3=0 


where  wo  =  wn  =  0.5  and  wj  =  1  for  j  ^  0,  N .  Now,  setting  r  =  n  for  i  =  0, 1, . . . ,  JV,  we  get 
the  non-linear  eigenvalue  problem 


{I-M{lm)\A  =  0,  m  —  0, 1, 2, ... ,  (7) 

where  M(7m)  is  a  square  matrix  of  size  2(JV+l)x2(iV  +  l)  and  A  =  [Ar  A (j)]T  is  the  column 
vector  consisting  of  the  unknowns  Ar<rn{ri)  and  A^m(ri).  For  each  given  m,  we  denote 

rm  —  |0  <  7m  —  ^0  sj^max  ~~  ^  = 

as  the  set  of  allowable  solutions  where  nmax  =  maxo<r<a  n(r).  The  problem  (7)  has  solution  if 
and  only  if 

det[/-Af(7m)]  =  0.  (8) 


3.2  Root  finding  problem 

The  solutions  of  the  determinant  equation  (8)  can  be  found  by  either  Newton-Raphson’s  method 
or  Ridder’s  method.  These  methods  are  both  implemented  for  this  report.  Let  go  from  0  to 
T  =  k0 i/max  n2  —  n2.  When  both  the  real  part  and  imaginary  of  the  determinant  change  their 
sign  in  one  step,  we  use  the  Newton-Raphson  method  or  Ridder’s  method  to  search  for  the  zero 
of  the  determinant  in  this  interval.  In  the  Newton-Raphson,  the  derivative  of  M  can  be  expressed 
in  a  closed  form  and  we  have 


~(*+l)  =  ryW  _ 
Jm  Im 


trace  j  [M{  7$ )]  ~ 1  (7m) )  } 


k  —  0, 1, 2, , 


(9) 
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where  trace  is  the  trace  of  a  matrix,  and  7m'  is  the  midpoint  of  the  interval  in  consideration.  In 
(9)  we  need  to  compute  the  matrix 

dM 


aim 


One  can  efficiently  do  this  by  solving 


m(^)b(^)  = 

ajm 


(10) 


To  simplify  the  notations,  we  denote  the  above  problem  as  MB  =  D.  We  use  the  following 
procedure  [6] 

Compute  PM  =  LU. 
for  k  =  1  :  m 

Solve  Ly  =  PD 
Solve  Ubk  =  y 

end 

where  Dk  is  the  fcth  column  of  D . 

4  Numerical  experiments 

In  this  section,  we  present  several  numerical  examples  which  illustrate  the  convergence  and  ac¬ 
curacy  of  the  integral  equation  formulation.  We  first  consider  the  step-index  fiber  since  its  exact 
solutions  are  available.  It  can  be  seen  from  the  computed  data  that  the  integral  formulation  is 
very  accurate  with  even  the  simplest  discretization  method  (trapezoidal  rule).  In  the  following 
numerical  examples,  the  parameters  are  unitless.  We  choose  operating  wavelengths  in  terms  of 
the  radius  of  a  fiber.  In  each  example,  we  plot  the  absolute  value  of  the  normalized  field  solution 
of  Er  and  of  Er  as  a  function  of  r  and  as  a  function  of  All  results  are  computed  with 

Matlab. 


4.1  Step- index  fiber 

Consider  a  step- index  fiber  with 

no  =  1.6,  nc  —  1.0,  a  =  1,  A  =  2a. 
In  this  case,  we  can  obtain  the  transcendental  equation 


(Jin 

\uJrr 


(«)  ,  K'mM  \  ( n0  JinM  ■ 


(u)  +  wKm(w)J  \n2uJm{u)  '  wKm(w) 


0( 


+ 


where 


u  —  a\J  /cqTIq  —  /?2,  w  =  a\J fi2  —  k^nl,  ko  — 


(11) 


-)  +  +  (12) 

))  \U2  W2 )  \TIq  u2  w2  ) 


In  Tables  1,  2,  and  3,  we  compare  the  solutions  7,,,  to  (8)  with  the  “exact”  solutions  to  the 
transcendental  equation  (12).  In  each  table,  N  denotes  the  number  subintervals  in  [0, a].  If 
m  —  0,  we  obtain  7 te  and  7 tm-  In  Figures  2  and  3,  we  plot  the  electric  field  distributions 
and  their  intensity  distributions,  respectively.  It  is  evident  from  the  tables  and  figures  that 
solutions  to  integral  equations  are  very  accurate.  The  convergence  rate  can  be  improved  when 
more  sophisticated  quadratures  are  applied. 
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N 

7TM 

7TE 

32 

2.18506665 

2.55757995 

64 

2.18337629 

2.55601272 

128 

2.18295441 

2.55562134 

512 

2.18282263 

2.55549906 

1024 

2.18281604 

2.55549294 

exact 

2.18281384 

2.55549091 

Table  1:  Exact  solutions  and  integral  equation  solutions  for  the  TE  and  TM  modes  (mode  number 
m  =  0). 


Exact  and  numerical  solutions  for  TE  and  TM 


Figure  2:  Exact  and  Nystrom  solutions  for  TE  ( Er )  and  TM  ( E cases.  The  number  of  subin¬ 
tervals  used  in  the  Nystrom  method  is  N  =  128. 


Figure  3:  Intensity:  \Er\ 2  and  |i^|2. 


4  NUMERICAL  EXPERIMENTS 


9 


N 

71,2 

71,i 

32 

0.57962508 

3.34572448 

64 

0.57603586 

3.34467642 

128 

0.57513626 

3.34441567 

512 

0.57491122 

3.34435058 

1024 

0.57485495 

3.34433431 

exact 

0.57483619 

3.34432888 

Table  2:  Exact  and  integral  equation  7  for  the  hybrid  modes  (mode  number  m  =  1). 


N 

72,1 

32 

2.20503141 

64 

2.20316831 

128 

2.20270322 

512 

2.20258699 

1024 

2.20255793 

exact 

2.20254824 

Table  3:  Exact  and  integral  equation  7  for  the  hybrid  modes  (mode  number  m  =  2). 

4.2  Multilayered  fiber 

Consider  the  following  multilayered  dielectric  fiber  immersed  in  air  with  the  air  core  and  12  layers 
of  dielectric  pairs,  which  can  schematically  be  written  as  n™*e(riHnL) I2nlfr’  The  fiber  parameters 
are 

nair  =  1,  nH  =  4.6,  nL  =  1.5,  (  , 

dcore  =  1,  dH  -  0.25,  dL  =  0.75.  1  ' 

Here,  =  1  denotes  the  air-core  region  and  =  1  denotes  the  air  region  exterior  to  the 
fiber.  The  fiber  radius  is  dmax  =  1  +  12(.25  +  .75)  —  13.  The  operating  wavelength  is  A  =  15.5. 
Hence,  if  the  unit  of  the  radius  is  in  millimeters  (mm),  the  corresponding  wavelength  is  15.5  mm 
which  is  in  the  microwave  region.  In  Figures  4  and  5,  we  plot  the  solutions  Er  and  E $  and  their 
intensity  distributions  for  m  =  1  and  7  =  1.03184.  Here  we  use  the  following  grid  sizes 

hcore  =  dcore/10,  hfj  —  dff/15,  hi,  =  djr//20. 

We  note  that  most  of  the  field  is  concentrated  in  the  air  core  and  the  first  layer  (0  <  r  <  2).  The 
ability  to  guide  waves  in  the  air  core  is  of  importance  since  non-linear  effects  and  propagation  loss 
are  reduced  significantly. 

4.3  Quarter-wavelength  fibers 

We  now  consider  quarter-wavelength  multilayered  fibers,  that  is,  their  radii  are  a  quarter  of  the 
operating  wavelength.  This  has  no  relation  to  the  quarter- wave  dielectric  mirrors  or  stacks  [5]. 
The  thickness  of  the  high-index  layer  and  the  low-index  layer  can  be  identical  or  different.  In 
the  following  examples,  we  solve  for  propagating  wave  whose  mode  number  is  m  =  1.  It  can  be 
seen  from  the  plots  of  intensity  distributions  that  most  of  propagating  field  is  concentrated  in 
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the  air  core  and  its  neighboring  layers.  This  suggests  that  low  power-loss  is  possible  for  these 
fibers.  In  the  microwave  region,  these  fibers  have  radii  between  7  mm  and  19  mm.  Furthermore, 
by  comparing  the  corresponding  solutions  gamma  (m  =  1)  in  each  case,  for  example,  Figures6 
and  ??,  8  and  22,  12  and  26  fibers  with  different  layer  thickness  can  confine  wave  in  their  air  core 
more  efficiently  than  those  with  same  layer  thickness. 

4.3.1  Different  layer  thickness 

Consider  fibers  with  6,  12,  and  18  layers  of  dielectric  pairs,  that  is,  na?re(ntfn£)fcnI?r  where 
k  =  6, 12, 18.  Fiber  parameters  are 

nair  =  1,  riH  -  4.6,  nL  =  1.5, 
dcore  =  1)  dfj  —  0.25,  d,L  =  0.75. 

In  Table  4,  solutions  7  to  (7)  for  m  =  1  are  listed  for  quarter-wavelength  multilayered  fibers  whose 
high-index  and  low-index  thicknesses  are  different. 


Number  of  layers 

7 

6  layers 

A  =  4  [dCOre  +  6  (d#  +  d£)]  —  28 

0.54353618 

12  layers 

A  =  4[dcore  +  12(dtf  4-  djr,)]  =  52 

0.03063851 

0.29596912 

18  layers 

A  =  4  [dCore  +  18(d#  +  d/,)]  =  76 

0.03105617 

0.20275634 

Table  4:  Solutions  7  for  multilayered  quarter- wavelength  fibers  with  different  number  of  layers. 
High-index  and  low-index  thicknesses  are  different. 
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Figure  6:  Normalized  intensity  profile  in  the  6-layer  fiber  of  radius  7  mm  with  different  refractive- 
index  thicknesses:  7  —  0.54353618. 


E,  n  =  1 


E,n  =  1 


Figure  7:  Normalized  intensity  distribution  in  the  6-layer  fiber  of  radius  7  mm  with  different 
refractive-index  thicknesses:  7  =  0.54353618 
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Figure  8:  Normalized  intensity  profile  in  the  12-layer  fiber  of  radius  13  mm  with  different 
refractive-index  thicknesses:  7  =  0.03063851 


Er,n  =  1 


Efn  =  1 


Figure  9:  Normalized  intensity  distribution  in  the  12-layer  fiber  of  radius  13  mm  with  different 
refractive-index  thicknesses:  7  =  0.03063851 
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Mode  number  n=1 


Figure  10:  Normalized  intensity  profile  in  the  12-layer  fiber  of  radius  13  mm  with  different 
refractive-index  thicknesses:  7  =  0.29596912 


Figure  11:  Normalized  intensity  distribution  in  the  12-layer  fiber  of  radius  13  mm  with  different 
refractive-index  thicknesses:  7  =  0.29596912 
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Mode  number  n=1 


Figure  12:  Normalized  intensity  profile  in  the  18-layer  fiber  of  radius  19  mm  with  different 
refractive-index  thicknesses:  7  =  0.031056172 


E,n: 


E,n  =  1 


Figure  13:  Normalized  intensity  distribution  in  the  18-layer  fiber  of  radius  19  mm  with  different 
refractive-index  thicknesses:  7  =  0.03105617 
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Figure  14:  Normalized  intensity  profile  in  the  18-layer  fiber  of  radius  19  mm  with  different 
refractive-index  thicknesses:  7  =  0.20275634 


Figure  15:  Normalized  intensity  distribution  in  the  18-layer  fiber  of  radius  19  mm  with  different 
refractive-index  thicknesses:  7  =  0.20275634 
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4.3.2  Same  layer  thickness 

Consider  fibers  with  6,  12,  and  18  layers  of  dielectric  pairs,  that  is,  ^a?re(n^n^)fcnafr  where 
k  =  6, 12, 18.  We  assume  that  the  high-index  and  low-index  layers  have  the  same  thickness.  Fiber 
parameters  are 


nair  =  1,  nH  =  4.6,  nL  =  1.5, 
dcore  “  1?  dff  —  0.5,  d[j  “  0.5. 

Solutions  7  with  m  =  1  for  different  fibers  are  listed  in  Table  5. 


Number  of  layers 

7 

6  layers 

A  =  4[<fCOre  +  6(dff  +  di)}  =  28 

0.26877105 

0.34457807 

0.71278887 

12  layers 

A  =  4[dcore  +  12  (<1h  +  <1l)\  =  52 

0.21106533 

0.38696444 

18  layers 

A  =  4[c?core  +  18  (dn  +  di)\  =  76 

0.14845343 

0.26511726 

Table  5:  Solutions  7  for  multilayered  quarter- wavelength  fibers  with  different  number  of  layers. 
High-index  and  low-index  thicknesses  are  equal. 


Mode  number  n=1 


Figure  16:  Normalized  intensity  profile  in  the  6-layer  fiber  of  radius  7  mm  with  same  refractive- 
index  thicknesses:  7  =  0.26877105 
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Figure  17:  Normalized  intensity  distribution  in  the  6-layer  fiber  of  radius  7  mm  with  same 
refractive-index  thicknesses:  7  =  0.26877105 


Mode  number  n=1 


Figure  18:  Normalized  intensity  profile  in  the  6-layer  fiber  of  radius  7  mm  with  same  refractive- 
index  thicknesses:  7  =  0.34457807 
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Er,n  =  1 


Ef"  =  1 


Figure  19:  Normalized  intensity  distribution  in  the  6-layer  fiber  of  radius  7  mm  with  same 
refractive-index  thicknesses:  7  =  0.34457807 


Mode  number  n=1 


Figure  20:  Normalized  intensity  profile  in  the  6-layer  fiber  of  radius  7  mm:  7  =  0.71278887 


4  NUMERICAL  EXPERIMENTS 


20 


Ef,  n  =  1 


*vn  =  i 


Figure  21:  Normalized  intensity  distribution  in  the  6-layer  fiber  of  radius  7  mm  with  same 
refractive-index  thicknesses:  7  =  0.71278887 


Mode  number  n=1 


Figure  22:  Normalized  intensity  profile  in  the  12-layer  fiber  of  radius  13  mm:  7  —  0.21106533 
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Er,  n  =  1 


E,,  n  =  1 

T 


Figure  23:  Normalized  intensity  distribution  in  the  12-layer  fiber  of  radius  13  mm  with  same 
refractive-index  thicknesses:  7  =  0.21106533 


r 


Figure  24:  Normalized  intensity  profile  in  the  12-layer  fiber  of  radius  13  mm  with  same  refractive- 
index  thicknesses:  7  =  0.38696444 
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Figure  25:  Normalized  intensity  distribution  in  the  12-layer  fiber  of  radius  13  mm  with  same 
refractive-index  thicknesses:  7  =  0.38696444 


Mode  number  n=1 


Figure  26:  Normalized  intensity  profile  in  the  18-layer  fiber  of  radius  19  mm  with  same  refractive- 
index  thicknesses:  7  =  0.14845343 
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Figure  27:  Normalized  intensity  distribution  in  the  18-layer  fiber  of  radius  19  mm  with  same 
refractive-index  thicknesses:  7  =  0.14845343 


Mode  number  n=1 


Figure  28:  Normalized  intensity  profile  in  the  18-layer  fiber  of  radius  19  mm  with  same  refractive- 
index  thicknesses:  7  =  0.26511726 
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Figure  29:  Normalized  intensity  distribution  in  the  18-layer  fiber  of  radius  19  mm  with  same 
refractive-index  thicknesses:  7  =  0.26511726 

5  Summary  and  plans 

In  this  final  report,  we  have  summarized  the  progress  made  during  the  period  September  1,  2002- 
September  30,  2003.  We  derive  the  integral  equations  for  electric  fields  propagating  in  a  fiber  of 
arbitrary  cross  section  and  arbitrary  refractive  index. 

Mathematical  analysis  of  conventional  optical  fibers  is  rich  and  widely  available  in  literature, 
for  example  [3,  9,  11]  and  references  therein.  It  helps  the  advancement  of  optical  fiber  industry 
in  understanding  and  designing  fibers  for  data  transport  and  telecommunications.  In  contrary, 
to  our  best  knowledge,  there  is  little  of  rigorous  mathematical  study  of  Bragg  fibers  due  to  the 
fabrication  difficulty.  Recently,  the  publications  of  A  dielectric  omnidirectional  reflector  [4],  An 
all-dielectric  coaxial  waveguide  [8],  and  External  reflection  from  omnidirectional  dielectric  mirror 
fibers  [7]  in  Science  report  successful  fabrications  of  omnidirectional  reflectors  and  multilayered 
fibers  for  optical  wavelengths.  These  articles  create  a  new  interest  in  studying  low-loss  Bragg 
fibers  at  microwave  and  millimeter  wavelengths  for  radar  applications.  Erom  previous  examples 
in  Section  4,  multilayered  fibers  of  radii  between  7  mm  and  19  mm  concentrate  wave  propagation 
within  their  air  core  for  wavelengths  in  the  microwave  range.  Hence,  they  can  be  very  useful  in 
radar  technology. 

A  thorough  analysis  of  the  integral  operators  described  in  Section  2  is  needed.  Their  spectral 
properties  can  be  studied  with  functional  analysis.  A  priori  knowledge  of  the  band  structure  of 
the  integral  operators  will  be  very  important  in  designing  accurate  and  fast  numerical  algorithms. 
Possible  problems  such  as  the  existence  of  spurious  modes  can  also  be  addressed  and  remedied  if 
one  can  predict  the  number  and  location  of  propagating  modes.  Error  and  convergence  analysis 
of  the  presented  numerical  method  needs  to  be  done  rigorously.  This  study  is  very  important, 
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especially  when  the  method  is  applied  in  optimal  designs.  We  plan  to  perform  these  analyses 
rigorously. 

For  radar  applications,  propagating  power  of  a  fiber  is  probably  the  most  important  property 
that  has  great  impact  on  the  fiber  usefulness.  Hence,  optimizing  its  power  over  parameters  such  as 
number  of  layers,  layer  thickness,  and  refractive  index,  is  of  great  interest.  We  plan  to  determine 
an  appropriate  objective  function  to  be  optimized  by  an  accurate  and  fast  optimizing  algorithm 
which  can  be  parallelized  for  different  computer  platforms. 

Acknowledgement 

The  first  author  would  like  to  thank  Dr.  William  Kent  and  Dr.  Gregory  Wilson  at  Mission 
Research  Corporation  in  Dayton  for  helpful  discussions  and  comments.  Dr.  Kent  brought  to  my 
attention  the  valuable  applications  of  low  loss  fibers  in  the  radar  world. 


References 

[1]  F.  Abeles,  Investigations  on  the  propagation  of  sinusoidal  electromagnetic  waves  in  stratified 
media .  Application  to  thin  films,  Ann.  Phys.  5  (1950),  596-. 

[2]  _ ,  Investigations  on  the  propagation  of  sinusoidal  electromagnetic  waves  in  stratified 

media.  Application  to  thin  films.  II  Thin  films,  Ann.  Phys.  5  (1950),  706-. 

[3]  A.  Bamberger  and  A.S.  Bonnet,  Mathematical  analysis  of  the  guided  modes  of  an  optical 
fiber ,  SIAM  J.  Math.  Analys.  21  (1990),  no.  6,  1487-1510. 

[4]  Yoel  Fink,  Joshua  Winn,  Shanhui  Fan,  Chiping  Chen,  Jurgen  Michel,  John  Joannopoulos, 
and  Edwin  Thomas,  A  dielectric  omnidirectional  reflector ,  Science  282  (1998),  1679-1682. 

[5]  Sh.A.  Furman  and  A.V.  Tikhonravov,  Basics  of  Optics  of  Multilayer  Systems ,  Basics  of, 
Editions  Frontieres,  Singapore,  1992. 

[6]  Gene  H.  Golub  and  Charles  F.  Van  Loan,  Matrix  Computations,  third  ed.,  Johns  Hopkins 
University  Press,  Baltimore,  1996. 

[7]  S.D.  Hart,  G.R.  Maskaly,  B.  Temelkuran,  P.H.  Prideaus,  J.D.  Joannopoulos,  and  Y.  Fink, 
External  reflection  from  omnidirectional  dielectric  mirror  fibers,  Science  296  (2002),  510-513. 

[8]  M.  Ibanescu,  Y.  Fink,  E.L.  Thomas,  and  J.D.  Joannopoulos,  An  all- dielectric  coaxial  waveg¬ 
uide,  Science  289  (2000),  415-419. 

[9]  P.  Joly  and  C.  Poirier,  A  numerical  method  for  the  computation  of  electromagnetic  modes  in 
optical  fibers,  Math.  Meth.  in  Appl.  Sciences  22  (1999),  389-447. 

[10]  Rainer  Kress,  Linear  Integral  Equations ,  second  ed.,  Springer- Verlag,  New  York,  1999. 

[11]  Rolando  Magnanini  and  Fadil  Santosa,  Wave  propagation  in  a  2-D  optical  waveguide ,  SIAM 
J.  Appl.  Math.  61  (2000),  no.  4,  1237-1252. 


REFERENCES 


26 


[12]  G.  Ouyang,  Y.  Xu,  and  A.  Yariv,  Theoretical  study  on  dispersion  compensation  in  air-core 
Bragg  fibers,  Optics  Express  10  (2002),  no.  17,  899-908. 

[13]  Y.  Xu,  G.  Ouyang,  R.  Lee,  and  A.  Yariv,  Asymptotic  matrix  theory  of  Bragg  fibers ,  J. 
Lightwave  Tech.  20  (2002),  no.  3,  428-440. 

[14]  A.  Yariv  and  P.  Yeh,  Optical  Waves  in  Crystals.  Propagation  and  Control  of  Laser  Radiation, 
John  Wiley  &  Sons,  New  York,  1984. 

[15]  P.  Yeh,  A.  Yariv,  and  E.  Marom,  Theory  of  Bragg  fiber,  J.  Opt.  Soc.  Am.  68  (1978),  1196- 

1201. 


